CU-TP-1187 

Interactions of massless monopole clouds 

Christopher M. Miller* and Erick J. Weinberg''" 
Physics Department, Columbia University, New York, New York 10027 



Abstract 

Some spontaneously broken gauge theories with unbroken non-Abelian gauge groups contain 
massless magnetic monopoles that are realized classically as clouds of non-Abelian field surrounding 
one or more massive monopoles. We use moduli space methods to investigate the properties of 
these massless monopole clouds. We show that the natural metric on the Nahm data for a class of 
SU(2M + 2) solutions with 2M massive and M{2M — 1) massless monopoles can be obtained from 
that for a simpler class of SU(M + 1) solutions. For the M = 1 case, we show that the Nahm data 
metric is isomorphic to the metric for the moduli of the BPS solutions, thus verifying a previously 
conjectured result. For M = 2 we use our results and the moduli space approximation to obtain 
an effective Lagrangian for an axially symmetric class of solutions. Using this Lagrangian, we 
study the interactions between the two types of clouds that appear. We show that, although the 
static spacetime field configurations suggest that the clouds might be rather diffuse, in scattering 
processes they behave as if they were relatively thin hard shells. 
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I. INTRODUCTION 



= 4 supersymmetric Yang-Mills theory is believed to possess an electric-magnetic 
duality symmetry. When the gauge group is maximally broken, to an Abelian subgroup, 
this duality interchanges the massive electrically charged particles corresponding to the el- 
ementary fields with magnetically charged states arising from classical soliton solutions. 
If the unbroken symmetry has a non- Abelian component, then there exist massless parti- 
cles carrying electric-type charge. Although there are no isolated massless solitons, certain 
multisoliton solutions have degrees of freedom that are naturally interpreted as correspond- 
ing to massless magnetic monopoles. These are manifested as clouds of non- Abelian field 
that surround one or more massive monopoles and shield part of their non-Abelian mag- 
netic charge [1]. Our aim in this paper will be to explore the properties of these massless 
monopoles by studying the interactions between these clouds. 

To explain this in more detail, consider an SU(iV) gauge theory with an adjoint Higgs 
field whose asymptotic value can be brought into the form 



with Si < S2 < . . . Sat. If the Sj are all distinct, the gauge symmetry is broken maximally, to 
U(l)^~^, and there are N — 1 topological charges. If the asymptotic magnetic field is then 
written in the form — eijkVk Qm/t^^, with 



these topological charges are the n^. We will refer to an SU(A^) monopole solution with 
such charges as being an (ni,n2, . . . ,njv_i) solution. 

With this maximal symmetry breaking, one can identify A^— 1 fundamental monopoles [2], 
each carrying a single unit of one of the topological charges, with the mass of the kth being^ 
(47r/e)(sjfc+i — Sfc). Each of these has four degrees of freedom — three position variables and 
one U(l) phase. A BPS solution with arbitrary magnetic charges can be understood as being 
composed of appropriate numbers of the various species of fundamental monopoles and as 
living on a moduli space whose dimension is four times the total number of component 
monopoles. 

Our interest here is not in maximal gauge symmetry breaking, but rather in the alter- 
native possibility, where $ has degenerate eigenvalues. The unbroken symmetry is then 
enhanced to a non- Abelian group, and some of the fundamental monopoles become mass- 
less. It is instructive to follow the behavior of the classical solutions as this case is obtained 
from the maximally broken one by smoothly varying the eigenvalues of One finds that 
an isolated fundamental monopole solution goes over to the vacuum solution as it becomes 
massless. The behavior of multimonopole solutions is more complex [3]. If the total mag- 
netic charge remains purely Abelian, then the solution rapidly approaches its limiting form 
once the inverse of the smallest monopole mass becomes larger than the separations of the 
component monopoles. In addition, the moduli space of solutions and its naturally defined 
metric have smooth limits; in the examples where the metric for the case with nonmaximal 



For the remainder of this paper, we will assume that the gauge fields have been rescaled so as to set the 
gauge coupling e to unity. 



$ = diag(si,S2,...SAr) 




Qm = diag (ni, n2 - Ui, . . . , -hn-i) , 
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breaking has been found directly, it is indeed the hmit of the metric for the maximally 

broken case [1, 4]. 

If instead the magnetic charge has a non-Abelian component, none of these are the 
case. Not only do the moduli spaces not have smooth limits, but cases that are gauge 
equivalent in the massless limit have different dimensions for any finite monopole mass. 
These pathologies are clearly related to the long-range behavior of the non-Abelian fields 
and the non-normalizability of certain zero modes in the massless limit. A simple example 
of this arises in an SU(3) gauge theory. With maximal breaking, the unbroken symmetry is 
U(l) xU(l), and there are two species of fundamental monopoles. If two of the eigenvalues of 
$ are equal, the unbroken symmetry is SU(2) xU(l) and one of the fundamental monopoles is 
massless. In the maximally broken case, the moduli spaces of the (2,0) and the (2,2) solutions 
have dimensions four and twelve, respectively. Yet, when the breaking is nonmaximal an 
SU(2) gauge transformation can turn a (2,[0]) solution into a (2, [2]) solution, where the 
square brackets denote the massless species. On the other hand, the (2,[1]) solutions, which 
have a purely Abelian long-range field have a well-defined moduli space, with a metric that 
is a smooth limit of the (2,1) metric. 

In this paper we will restrict ourselves to configuration with purely Abelian magnetic 
charge. Several such solutions with a single massless monopole are known. These include 
solutions [5] with one massive and one massless monopole for S0(5) broken to SU(2) xU(l), 
as well as SU(4) (1, [!],!) solutions [6] and the SU(3) (2,[1]) solutions [7] referred to above, 
each with one massless and two massive monopoles. These SU(3) solutions will be of par- 
ticular importance in our analysis. Their properties were studied in considerable detail by 
Dancer and collaborators [7-10], and we will refer to them, and their (A'^ — 1, [A^ — 2], . . . , [1]) 
generalizations for SU(A^), as Dancer solutions. 

In all of these solutions there is a single non-Abelian cloud (spherical in the first case, 
ellipsoidal in the latter two) that encloses the massive monopole(s). Well inside this cloud, 
the magnetic field approximates the field, with both Abelian and non-Abelian magnetic 
components, that would be expected to arise just from the massive monopoles. The cloud 
effectively shields the non-Abelian components, so that outside the cloud one sees the purely 
Abelian field corresponding to the sum of the massive and massless monopole charges. The 
size of the cloud is determined by a single parameter, whose value has no effect on the total 
energy of the solution. 

Since our goal in this paper is to investigate the interaction between massless monopole 
clouds, we need solutions that have more than one cloud. One might have expected that 
these could be obtained simply by having more than one massless monopole. This turns 
out not to be so. For example, the SU(A") (1,[1],. . . ,[1],1) solutions contain — 3 massless 
monopoles, but only a single cloud, no matter how large A^ becomes [1]. [In fact, for A^ > 4 
these solutions are essentially embeddings of (1,[1],1) SU(4) solutions into the larger group.] 

A set of solutions that do have multiple clouds, and which we will focus on, are the 
(2, [2], [2], [2], 2) solutions — with four massive and six massless monopoles — in the theory 
with SU(6) broken to U(l) xSU(4) xU(l). The structure of these solutions was studied 
in Ref. [11]. They can be viewed as containing two SU(3) Dancer solutions, each with a 
"Dancer cloud" enclosing two massive monopoles, embedded in disjoint subgroups of the 
SU(6). In addition, there are two larger "SU(4) clouds", enclosing both Dancer clouds, that 
are somewhat analogous to the cloud in the SU(4) (1,[1],1) solutions. As we will describe 
in more detail later, these four clouds are characterized not only by individual cloud-size 
parameters, but also by a number of additional parameters specifying their orientations with 
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respect to the unbroken gauge group. Note that, even in this case, the number of clouds is 
less than the number of massless monopoles. 

A useful tool for studying the low-energy dynamics of multimonopole systems such as 
these is the moduli space approximation [12], which reduces the full field theory dynamics 
to that of a finite number of collective coordinates z"'. The latter are governed by the 
Lagrangian 

L^)^g^z''z\ (1.3) 

where Qah is the metric on the moduh space of BPS solutions. This method was used to 
study the (1,[1],1) SU(4) solutions [13]. Because the two massive monopoles lie in mutually 
commuting subgroups of the SU(4), one would not expect them to interact directly with each 
other. This turns out to be the case; indeed, they can pass though each other undeflected. 
On the other hand, they do interact with, and exchange energy with, the cloud. In these 
interactions the cloud acts as if it were a thin shell, with the monopole-cloud interactions 
concentrated in the times when the massive monopole positions coincide with the shell. 
At large times the massive monopoles and the cloud decouple from each other and evolve 
independently. The interactions between the cloud and the massless monopoles are similar 
in the SU(3) Dancer solutions [8, 9], which have the additional feature that the massive 
monopoles, being both of the same type, also interact directly with each other. 

A word of caution is in order here. The essential idea underlying the moduli space approx- 
imation is that for a slowly moving soliton fluctuations off of the moduli space are energet- 
ically suppressed. If the theory contains massless particles, this needs further examination, 
since excitation of these modes by radiation of massless particles is always energetically pos- 
sible. However, since the source of the radiation is proportional to the time derivative of the 
bosonic flelds, the radiation rate is expected to be small for low soliton velocities. This has 
been shown rigorously for configurations involving pairs of monopoles in theories where the 
massless gauge fields are all Abelian [14, 15]. Although the validity of the approximation 
has not been demonstrated rigorously for the case where there are massless non-Abelian 
gauge fields, explicit comparison of the predictions of the moduli space approximation with 
numerical evolution of the full field equations in a spherically symmetric example [16] show 
that they agree as long as the configurations are slowly varying. 

Of course, this method requires that moduli space metric be known. It can be obtained 
directly from the BPS solutions, if they are known explicitly. In some other cases indirect 
methods based on the mathematical properties of the moduli space can be used to obtain 
Qah [17-19]. However, for other cases it turns out to be easiest to resort to an alternative 
approach. The Atiyah-Drinfeld-Hitchin-Manin-Nahm (ADHMN) construction [20-23] is a 
powerful tool for obtaining BPS solutions. It is based on an equivalence between the Bogo- 
molny equation for the fields in three-dimensional space and an ordinary differential equation 
for a set of matrix functions of a single variable, known as the Nahm data. The moduli space 
of Nahm data has its own naturally defined metric. It has been shown for both the SU(2) 
theory [24] and for the case of SU(A'" -|- 1) broken to U(iV) [25], and is believed to be true 
in general, that the moduli spaces of Nahm data and of BPS solutions are isometric.^ [In 
particular. Dancer used this equivalence in his investigation of the dynamics of the SU(3) 
solutions, and worked with the Nahm data metric] In this paper we will assume that the 



^ In fact, we will demonstrate this equivalence for yet another example, the (1, [!],!) SU(4) solutions, in 
Sec. III. 
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equivalence of the two metrics holds as well for our SU(6) example, and will work with the 
Nahm data metric. 

The remainder of this paper is organized as follows. In Sec. II we review the relevant parts 
of the ADHMN construction, as well as the metric on the moduli space of Nahm data. We 
then show that if the metric for the SU(M+ 1) Dancer sohitions is known, then that for the 
(M, [M], [M],M) solutions of SU(2M + 2) broken to U(l) xSU(2M) xU(l) can be easily 
obtained. In Sec. Ill we illustrate this method by obtaining the metric for the (1, [!],!) SU(4) 
solutions and verify that the Nahm data metric thus obtained is in fact equivalent to the 
metric on the moduli space of BPS solutions, which had been found previously [1, 13, 26]. 
Next, in Sec. IV, we turn to the moduli space metric for the SU(6) (2, [2], [2], [2], 2) solutions. 
These are described by a total of 40 collective coordinates. Although the full metric can be 
obtained by the methods of Sec. II, the result would be rather unwieldy for exploring the 
nature of the cloud dynamics. We therefore consider a restricted problem, that of a lower- 
dimensional submanifold of axially symmetric solutions, for which the metric has a relatively 
simple closed form expression but which still has enough structure to allow us to investigate 
nontrivial cloud dynamics. In Sec. V we use the metric obtained in the previous section to 
explore the cloud-cloud interactions. Section VI contains some concluding remarks. There 
is an Appendix, which contains some of the details of the moduli space metric calculation. 

II. THE ADHMN CONSTRUCTION AND THE METRIC FOR 
(M, [M], . . . , [M], M) SOLUTIONS 

As was discussed in Sec. I, the moduli spaces of BPS solutions and of Nahm data are 
believed to be isometric. We work in this paper with the Nahm data metric, which is the 
more accessible of the two for the examples that we study. In this section we first review 

the essential elements of the ADHMN construction [20-23], emphasizing the points that 
arc relevant for the {M. [M]. . . . , M) solutions of SU(2A'/ + 2). We then obtain a general 
expression for the metric of these SU(2M + 2) solutions, and briefly discuss an asymptotic 
special case. 

A. Nahm data 

The basic elements in the ADHMN construction^ are the Nahm data, a quadruple of 
matrix functions T^(s) (/^ = 0, 1, 2, 3) that obey the Nahm equation, 

dT- i 

= ^ + ^[To, Ti] + -e,,k[Tj, T^] , t, j, k = 1,2,3. (2.1) 

For charge k solutions in an SU(2) theory with Higgs vacuum expectation value v, the T^{s) 
are k x k Hermitian matrices defined for —v/2 < s < v/2. For a general SU(A^) theory, the 
eigenvalues of the Higgs vacuum expectation value divide the range of s into A^ — 1 intervals, 
on each of which Eq. (2.1) must hold. (The boundary conditions at the ends of these intervals 
are somewhat involved; we describe them below for the cases we need.) The dimension of the 



^ For a fuller description of the ADHMN construction, including a discussion of how the spacetime fields 
are obtained from the Nahm data, see Ref. [27] . 
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varies from interval to interval, being determined on each by the corresponding magnetic 
charge. If the are of the same size, k x k, on two adjacent intervals, then there are 
additional "jumping data", forming a 2/c-component complex vector, associated with the 
boundary between these intervals. 

For later reference, note that if the T'^(s) satisfy the Nahm equation, then so do the 
defined by 

fi{s) = R,,Tj{s) + Bi, 

fo(s) = To{s) , (2.2) 

where Rij is an orthogonal matrix with determinant one. The T^(s) lead to a spacetime 
solution that is obtained from the original one by a combination of the spatial rotation 
specified by R and a spatial translation by the vector B. 

For SU(2M + 2) broken to U(l) xSU(2M) xU(l) the eigenvalues of the Higgs vacuum 
expectation value are Sl < Sq < Sr, with Sq being 2M-fold degenerate. These divide the 
range of s into a "left" interval [sl,so], a "right" interval [so,Si?], and 2M — 1 intervals of 
zero width ai s — sq. These correspond to two species of massive monopoles, with masses 

Ml = Mso - sl) , Mr = A7r{sR - so) , (2.3) 

and 2M — 1 species of massless fundamental monopoles. For the (M, [M], . . . , [M], M) so- 
lutions we need two sets of M x M matrices, Tl^{s) and T^{s), defined on the left and right 
intervals, respectively.^ These obey Eq. (2.1) subject to the boundary condition that, for 
M > 1, the Tj^ (Tj^) have poles at (sr) with the residues forming an M-dimensional irre- 
ducible representation of SU(2). Except for these poles, the and must be everywhere 
nonsingular. 

Because the magnetic charge is the same on adjacent intervals, there are jump data 
associated with each of the 2M coincident boundaries at sq. We write the data associated 
with the Fth boundary as a 2M-component vector a^^, where r = 1, 2, . . . , M and a — 1,2. 
These jump data are required to satisfy the constraint 

- [Tn'o)i, - [Tfi^o)i, - I E Csi^^Uai . (2.4) 

If we assemble the jump data into a 2M x 2M matrix A with Aar,F — ^ar-' define an 
M X M matrix 

(T^U^l^Oar (2.5) 
^ Fa 

and a 2M x 2M matrix 



The spacetime fields are obtained from sums of integrals over the various intervals, with the integrands 
obtained by solving a differential equation involving the Nahm matrices on the corresponding interval. 
Because the integrals for the zero-width intervals vanish, the corresponding Nahm matrices have no effect 
on the spacetime fields. For more details, see Ref. [27]. 
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the jump data constraint becomes simply K = AA^ . For later reference, it is important to 
note that this implies that the eigenvalues of K must be positive. The general solution of 
this constraint is 

A = K^/^V, (2.7) 

where V is unitary. The freedom in the choice of V reflects the existence of an unbroken 
U(l)xSU(2M) subgroup of the original SU(2M + 2) gauge group.^ 

The crucial observation for us is that the dimensions and boundary conditions for the 
and are the same as those for the Nahm data of the SU(M + 1) Dancer solution. Apart 
from the gauge orientation angles and phases encoded in V, the new information specific 
to the SU(2M + 2) problem — i.e., the parameters that describe the SU(2M) clouds — 
enters only through T4, which can be chosen to be any Hermitian matrix, subject only to 
the constraint that the eigenvalues of K must be positive. Thus, if the Nahm equation has 
already been solved for the SU(M + 1) Dancer problem, the only additional data needed for 
the (M, [M], . . . , [M], M) solutions are the jump data, which are given by Eq. (2.7). 



B. Gauge action 

In addition to the spacetime symmetries described by Eq. (2.2), the Nahm equations have 

a set of invariances that are analogous to, although distinct from, the gauge transformations 
on the spacetime fields. If g{s) is a unitary matrix of appropriate dimension, the Nahm 
equation (2.1) and the jump equation (2.4) are invariant under^ 

T,{s) ^ f,{s)^g{s)T,{s)g-\s)+iS,o^g-\s), 

When considering spacetime gauge transformations one distinguishes between local gauge 
transformations, which approach the identity at spatial infinity, and global gauge transfor- 
mations, which are nontrivial as r — oo. While the former simply refiect the presence of re- 
dundant field components, the latter correspond to symmetries that are related to conserved 
gauge charges. When the gauge symmetry is unbroken, the number of global gauge trans- 
formations is equal to the dimension of the gauge group. If the symmetry is spontaneously 
broken, as it must be when monopoles are present, only those global gauge transformations 
that leave the asymptotic Higgs field invariant (i.e., those in the unbroken gauge group) lead 
to normalizable zero modes about a static solution, and only these correspond to physical 
motions on the moduli space. 

Similarly, we can distinguish between local gauge actions on the Nahm data, for which 
g{s) = / at both boundaries, and global gauge actions, which have g ^ I at one or both 
boundaries. Any of the latter that act on pole terms in the Nahm data will lead to non- 
normalizable modes, and therefore do not contribute to the moduli space dynamics. This 
means that for the (M, [M], . . . , M) solutions, which have poles at both boundaries, the 
only relevant global gauge actions are those that leave the pole terms invariant. These are 
proportional to the unit matrix and are of the form g{s) = e^^^^^I. Because an s-independent 



^ In the next subsection we will see how the effects of the remaining unbroken U(l) factor are manifested. 
^ Such transformations are often used to make To(s) vanish identically. 



7 



gauge action proportional to the unit matrix would have no effect on the Nahm data, it is 
sufficient to consider the case where x vanishes at one boundary, but not at the other; this 
leads to a zero mode corresponding to an unbroken U(l). The zero modes corresponding to 
the other unbroken generators do not arise from gauge actions, but instead correspond to 
variations of V . 

For the SU(M + 1) Dancer solutions there must be a pole at one boundary, but the 
only constraint at the other is that the Nahm data be nonsingular. There are then a total 
of independent normalizable global gauge zero modes, corresponding to the generators 
of the unbroken U(M). These can be obtained from gauge actions for which g is an M- 
dimensional unit matrix at the pole and proportional to one of the generators of U(M) at 
the other boundary. 

C. The moduli space metric 

The moduh space of Nahm data is the space of solutions of the Nahm and jump equa- 
tions, but with solutions that are related by local gauge actions considered equivalent. The 

coordinates on this space are the collective coordinates Za- Its tangent space at a given point 
is spanned by the variations of the Nahm data that preserve the Nahm and jump equations 
and that are orthogonal to the variations due to local gauge actions. 

For our solutions, the Nahm data consist of two quadruples of Nahm matrices, one for 
each interval, and jump data at the boundary at Sq. We can write these collectively as 
T = {T^(s), T^(s), A}. The tangent vectors to the Nahm data moduli space must have a 
similar structure, and so can be written SkS y = {Y'^(s), F}. The inner product of 

two such vectors is defined to be^ 

{y.y)= r '^Yi:{s)Yl^{s)+ r 1^Y^{s)yI\s) + \i^{YY"^ + Y'Y^). (2.9) 

An infinitesimal local gauge action is specified by a Hermitian matrix function A(s) that 
is everywhere continuous and vanishes at and sr. We denote its values on the left and 
right intervals by A^(s) and A^(s), and define A(so) = A°. Its action on the Nahm data 
defines a vector D^gauge with 



Note that in the trace in the last term the indices run over 2M values, whereas in the first two terms the 
traces are of M x M matrices. 
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(2.10) 



where Dj = i[Tj{s), ] and Dq — d/ds + i[To{s), ]. In order that y be orthogonal to 3^gauge 
for any choice of A, we must require that 







S7^ So, 



= 



yoiso)-Y,^{so)l^ + -{YA^-AY^), 



(2.11) 
(2.12) 



By analogy with the corresponding constraint on the variations of the spacetime fields, we 
will refer to these as the background gauge conditions. 

A basis for the tangent space is given by a set of vectors of the form 



an 



dA 



dz" 

- i[Aaiso)^A] 



(2.13) 



with the gauge action Aa{s) chosen so that is in background gauge. The metric on the 
moduli space is then defined to be^ 



ds' = gab dz'^dz" = in{ya, yb) dz'^dz' . 



o.jJ> 



(2.14) 



D. The metric for the (M, [M], . . . , M) solutions of SU(2M + 2) 

As we have just seen, to calculate the moduli space metric one must first vary the Nahm 
data with respect to the coordinates, and then find a gauge action A(s) that brings the 
resulting tangent vector into background gauge. The one part of this procedure that is not 
completely straightforward is the solution of the differential equation for A(s) that is implied 
by Eqs. (2.11) and (2.13). However, if this has already been done for the SU(M+ 1) Dancer 
solutions, the determination of the moduli space metric for the (M, [M] , . . . , M) solutions of 
SU(2M+2) reduces to an algebraic problem. To see this, first note that the Dancer problems 
give us two sets of basis vectors, Y^J" and 1^^^, that satisfy Eq. (2.11) on their respective 
domains. Each of these sets includes vectors corresponding to the global U(M) gauge 
freedom. These are of the form Yj^ = D^Xf ^iid Yj^ = D^xf (/ = 2, . . . M^), where Xf 
{xf) vanishes at sl {sr) and is nonzero and proportional to one of the U(M) generators at 
So- Because the Dancer basis vectors satisfy Eq. (2.11), 

D^D^Xf=D^D^Xf^O. (2.15) 

Now suppose that the coordinates for the (M, [M], . . . , M) solutions arc chosen so that 
one subset are those originating with the left Dancer problem, a second subset are those 



^ The factor of An here is chosen to make the normahzations of the Nahm data metric and the BPS solution 
metric the same; it can be easily checked by noting the coefficient of the terms quadratic in the center-of- 
mass position. 
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from the right Dancer problem, and the remainder are associated only with the jump data. 
The tangent vector corresponding to the coordinate z"" can then be written in the form 



V = I Y^^ + D \^ Y^^ + D \^ 



DR 
an 



dA 
dz"- 



z(A>/2)^ 



(2.16) 



where Yj^^ and Yj^^ are the vectors from the left and right Dancer problems.^ Although 
these Dancer vectors already satisfy the background gauge condition on their respective 
domains, an additional gauge action may be needed to satisfy Eq. (2.12) at sq. Its gauge 
function A must vanish at and and, in order to maintain the background gauge 
condition on the two intervals, it must satisfy D^DnK = for s ^ Sq. It is therefore a linear 
combination of Dancer global gauge modes, with 



A 
A 



L _ L L 



4xf 



[2.17) 



Thus, for each y^, there are a total of 2M^ constants to be determined. Requiring continuity 
of the gauge action at the boundary, A° = A^(so) = A^(so), gives algebraic equations. 
The background gauge condition at the boundary becomes 




dA 

dZa 



At 



(2.18) 



and gives more equations, thus determining the c^j and c^j, and hence y^. 

We can now evaluate the metric. Because the y^ satisfy the background gauge conditions 
of Eqs. (2.11) and (2.12) and, in addition, D^Y^^^^ = D^Y^^^^ = 0, many of the terms 
involving the gauge actions can be eliminated by integrations by parts. With the aid of 
Eq. (2.18), one eventually obtains 



g,, = 47r r ds Tr l^^^F^^ + 4n H ds Tr F^^^n^ + dlb 

JSL Jso 



DL , DR , 
iJab ^ iJab ^ {Jab ; 



(2.19) 



where g^^^ and g^j,^ are the metric components from the corresponding Dancer solutions and 

dAdA^ dAdA^\ 



r7° 
Uab 



47rTr 



(^o) - i;r(5o) A0 + 27rTr 



+27riTr 



OA 



-A\Xl®h)-{Kl®h)A^^ 



dz" dz'' dz^ dz"- 
dA^' 



= 27rTr 



dz 

dA dA^ dA dA^ 



\ dz"- dz^ dz^ dz" ^ 
47rTrT4(AX + A^A°) 



+ 47rIV DoA^iso) - DoA^iso) 



A? 



(2.20) 



^ Of course, with coordinates chosen as described above, at most one of these Dancer vectors will be nonzero 
for a given 0". 
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contains the entire contribution from the jump data.^° These two equations are the main 
result of this section. 

E. Large SU(2M) clouds 

Before focusing on the specific examples of M = 1 and M = 2, it is worth commenting 
briefly on the case where the length scales in Tj^^ and are smaller than all those entering 
T4 by a factor of e ^ 1. This corresponds to the situation where the SU(2M) clouds are large 
compared to both the Dancer clouds and the separations between the massive monopoles. 
In this case Eq. (2.7) can be written as 

A={tI'^ ®h)V + (2.21) 

where 5A^ which contains all of the information about the Dancer data, is suppressed by a 
factor of e. Therefore, if Za is one of the Dancer coordinates, dA/dza = 0(e). Furthermore, 
by noting the A°T4 terms on the left-hand side of Eq. (2.18), we see that the Aa corresponding 
to these coordinates arc also suppressed by a factor of e. It follows from these facts that 
(1) if either a or 6 refers to a Dancer coordinate, then g^i^ is suppressed relative to gj^^^ or 
gi^^ and (2) if a and b both refer to jump coordinates, then all of the dependence of g^j, on 
Dancer parameters is through subleading terms. Hence, to leading order in e the moduli 
space Lagrangian separates into two parts, one depending only on the Dancer parameters 
and one depending only on the jump parameters. In other words, large SU(2M) clouds are 
effectively decoupled from both the massive monopoles and the Dancer clouds. 

III. (1,[1],1) SOLUTIONS IN SU(4) 

We will first consider the (1,[1],1) solutions for a theory with SU(4) broken to 
U(l)xSU(2)xU(l). This will not only serve as an illustration of our method but, because 
the metric on the space of BPS solutions is already known [1, 13, 26], it will also provide 
one more example supporting the conjecture that the moduli spaces of Nahm data and of 
BPS solutions are isometric. 

The "Dancer" solutions in this case are just embeddings of the unit SU(2) monopole, each 
with a four-dimensional moduli space. The Nahm data are numbers rather than matrices 
and are given (with a standard choice of the gauge action) on the left interval by Tq = 0, 

= -Xf, where the are the coordinates of the monopole center. Differentiating these 
with respect to the Xj" gives three Dancer tangent vectors 

= (3.1) 

that satisfy the background gauge condition without needing any compensating gauge action. 
The fourth tangent vector corresponds to a U(l) phase, and so must be of the form 

J, X 

vDL 
_^U(1) 



Although the middle term in the final expression for g'^j^ appears not to be symmetric under interchange 
of a and b, it actually is. This can be shown by an integration by parts and making use of the facts that 
the Ac obey Df^D^Ac = and vanish at sl and sr. 



3 



/7c 



(3.2) 
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In order that this be in background gauge, we need that d^K/ds^ 
ization by requiring that X'^(so) = 1 and x^^^l) = 0, we find that 



47r(g - Sl) 

m7 



^U(l) 



An 



0. Fixing the normal- 



(3.3) 



The Nahm data and tangent vectors for the right Dancer data are completely analogous. 
They are obtained simply by replacing L by i? in the above equations, except for a sign 
change that yields 



X^(^) 



47r(s - sr) 



M, 



R 



VDR 
^U(l) 



-5 



An 



M, 



R 



(3.4) 



In the absence of jump data, these two sets of four vectors would give a moduli space metric 



dsl + dsl = Ml dXl + Mr dX^ + 



M 

Mr 



dxl + 



M, 



R 



R 



(3.5) 



In the standard fashion, we can rewrite the positions in terms of center-of-mass and 
relative positions Xqm and R. We can also replace xl and xr by a global U(l) phase and a 
relative U(l) phase. The former, given by ^ = Xl + XRj corresponds to a simultaneous phase 
rotation of the two monopoles, and is described by a tangent vector { 



VDR 



A7r5^o/{ML+MR){l, 1} = An/iML+MR) D^{Ml x"", -Mrx""}- Note that although 



yDL 



and 



VDR 



correspond to pure gauge actions on their separate intervals, the combined 



vector is not a gauge action because the corresponding left and right gauge functions are 
not equal at s = sq. 

The relative U(l) phase = {MlXr ~ MrXl) / {Ml + Mr) corresponds to an orthogonal 
combination of vectors and, in the context of just the Dancer data, is a pure gauge action. 



We can therefore choose the gauge so that 



and 



vdr 



both vanish and the Qatp are 



given completely by the jump data term g^^ 

We now have to consider the contributions from the jump data. We start by defining 
T4 = b; examination of the spacetime solutions shows that b measures the size of the non- 
Abehan cloud. We then have 



X = 6/2 + R • <T = UKqU' 



(3.6) 



where Kq = diag {b + R,b — R), and can write the general solution for the jump data as 

A = C/Xo/Ve^^ , (3.7) 

where W, like U, is an SU(2) matrix. Neither the center-of-mass position nor the global 
U(l) phase ^totai appears in A, so the tangent vectors corresponding to these variables have 
no jump component and are specified completely by the and Y^^ inherited from the 
Dancer problems. It is easily verified that these vectors are in background gauge and that 
they are orthogonal to the vectors for the relative coordinates. 



^ In the notation of Eq. (2.7), W — U ^V. We have written A in this form to facihtate comparison with 
the results in Ref. [13]. 



12 



The calculation of the remaining metric terms is simplified by noting that, because the 
spatial rotations represented by U and the global U(2) symmetry represented by W and 
ijj are isometries, we can calculate the metric from the tangent vectors at a point with 
U — W — I, t/j — 0. At this point the tangent vector for the intermonopole separation R 
gets a jump data contribution 

dA 



1 



that combines with the contributions from the Dancer vectors to give a tangent vector y^i 
that is already in background gauge. For the cloud size parameter b we have 



dA 
'db 



1 



diag 



1 



1 



2 \VbTR' Vb^J ■ 
Because b does not enter the Dancer solutions, yt has no l^f ^ or Yf^ contribution. 



(3.9) 



It, 



too, satisfies the background gauge conditions without the need for a compensating gauge 
action. 

To obtain the remaining tangent vectors, which correspond to rotations and U(2) trans- 
formations, we first write an infinitesimal variation of A (with Kq held fixed) as 



dA 



-1/2 



(3.10) 



where the daj and d^j are the invariant one-forms for the rotational SO (3) and gauge SU(2) 
symmetries. The dai and da2 terms combine with contributions, due to rotations of R, 
from the left and right intervals to give background gauge tangent vectors. The /9i and 
[32 vectors have no contributions from these intervals, but are also in background gauge. 
However, the and ip vectors do not satisfy Eq. (2.12) and so must be supplemented 
by compensating gauge actions. As explained in Sec. II D, these gauge actions must have 
gauge functions of the form = c^x^y = c^X^- Continuity of A at Sq implies that 



= Ca- Using Eq. (2.18), we then find that 



Cfh 



Rfj, 



bfi 



4tt + 2bii 



(3.11) 



4:71 + 2b^ 

where ^ = M^Mr/ {Ml + Mf{) is the reduced mass. 

With the Aa thus determined, we can use Eqs. (2.19) and (2.20) to show that the metric 
for the eight-dimensional relative moduli space [i.e., with the center-of-mass motion and 
overall U(l) phase factored out] is 



ds'' = 



27r6 



(b'^-R^) 



dR^ + 



27r6 



(62 - i?2) 



db^ 



AnR 
b^-R^ 



dbdR 



+ [liR^ + 27r6) {dal + dal) + 2T:b{d(3l + d(3l) + A-Ky/h^ - R^{dai d(3i + da2 d(32) 



27r + 6/i 



27r6- 



2nfiR^ 
{2n + bfi) 



df3l + 



Stt^R 
{271 + bfi) 



dip dPs 



(3.12) 



This agrees with the metric on the moduli space of BPS solutions that was previously 



We have not included an term because its effect can be absorbed by a redefinition of P3; tfiis term 
would correspond to the Euler angle that leaves R invariant. 
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obtained. This thus provides another example where the moduli spaces of Nahm data and 
of BPS solutions are isometric, lending further support to the conjecture that this is true in 
general. 

For later reference, we note that when the angular momenta and charges all vanish, the 
system reduces to one governed by the Lagrangian 

IV. (2, [2], [2], [2], 2) SOLUTIONS IN SU(6) 
A. Nahm data 

The Nahm data are 2x2 matrix functions T'^(s) and T^{s) on the left and right intervals, 
respectively, plus jump data that are obtained from T'^(so), T^{sq), and 

T,=ph + q-T (4.1) 

by using Eqs. (2.6) and (2.7). The parameters in T4 determine the properties of the two 
SU(4) clouds. Examination of the spacetime solutions [11] shows that, roughly speaking, 
p + q and p — q (with q = |q|) determine the sizes of the clouds, while the direction of q 
specifies an orientation in the unbroken SU(4). The T^(s) are themselves the Nahm data 
for an SU(3) (2,[1]) Dancer solution. By an appropriate gauge action one can set Tq — 
and then write [8] 

^^ = |E4//(«)^/ + ^fl2, (4.2) 

where Af^ is an orthogonal matrix and the three f^^ = UlTjU^^ are a rotated set of Pauli 
matrices. The //(s) obey 

f^/f/f (4.3) 

and its two cyclic permutations. If we adopt the convention that /f < /| < /|, they are 
given in terms of Jacobi elliptic functions by 



Dl cn^^ [Dl{s - Sl)] 
sn^^[DL{s - Sl)] 

Dl dn,^[DL{s - sl)] 
sn«^ [DLis -Sl)] 



The requirement that //(s) only have a pole at sl imposes the conditions < kl < 1 and 
< Dl < 2K[kl)/{so — sl), where K{k) is the complete elliptic integral of the first kind. 
The T^{s) are similar, but with Dl and kl replaced by Dr and kr. 



This is verified most easily by comparing with the form given in Ref. [13]. For the angular and phase parts 
of the metric our expression, written in terms of angular velocities, is related to the one in Eq. (3.18) of 
that paper, given in terms of angular momenta, by a Legendre transformation. 
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FIG. 1: A geodesically complete submanifold illustrating the SU(3) Dancer solutions. The long 
straight lines correspond to the axially symmetric hyperbolic solutions, with two widely separated 
monopoles, while the short straight lines correspond to the axially symmetry trigonometric solu- 
tions. The limiting curved boundaries, which are not part of the manifold, correspond to SU(2) 
two-monopole solutions. 

The left and right Nahm data each contain eleven parameters: three center-of-mass vari- 
ables Rj, the three Euler angles in that specify the spatial orientation, the three angles 
needed to define the fj, and the elliptic function parameters D and k. The significance of 
the latter two is clarified by referring to the plot in Fig. 1. The change of variables [9] 

X = (2 - «2)d2 ^ 

y = -a/3«2d2 (45) 

maps the allowed range of D and n onto the lower right sextant of the plot (including the 
straight boundaries, but excluding the curved outer boundary, which is geodesically infinitely 
far from any point in the interior). By adjoining five other copies (corresponding to the 
other possible orderings of the /?), one obtains a geodesically complete two-dimensional 
manifold. Points far out on the long arms of the figure correspond to solutions with two 
well-separated massive monopoles, with the distance between the two approximately equal 
to D. The straight fines down the centers of the arms, on which k — 1, correspond to 
minimal Dancer cloud size, while the limiting curve corresponds to embeddings of SU(2) 
two-monopole solutions that can be thought of as having infinite Dancer clouds. The central 
point, where D = Q and n is undefined, corresponds to a solution with coincident massive 
monopoles and a minimal size cloud. On the short straight lines emanating from this point 
K — Q. Points on these lines correspond to solutions with coincident massive monopoles and 
clouds varying from minimal to infinite size [10, 28]. 

For K equal to or 1, the elliptic functions reduce to trigonometric or hyperbolic functions, 
respectively [7]. Two of the fj are then equal and the spacetime solution has an axial 
symmetry. When D — 0, all three of the fj are equal and the solution is spherically 
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symmetric. A straight trajectory passing from a k = line through the central point and 
out along the opposite k — 1 line is a geodesic of the Dancer metric. 



B. Reduction to cylindrical symmetry with vanishing charges 

As noted above, the left and right sets of Dancer data each contain eleven parameters. 
In addition, there is an overall U(l) phase associated with each set. These, plus the four 
parameters from the elements of T4 given in Eq. (4.1) and the sixteen moduli arising from the 
U(4) matrix V would seem to give a total of 44 moduli. This cannot be correct, because a 
solution with ten monopoles should lie on a 40-dimensional moduli space. The discrepancy is 
resolved by noting that there is a U(2) subgroup of the U(4) whose effect is gauge equivalent 
to that obtained by simultaneously rotating the U(l) phases and SU(2) orientations of the 
two Dancer solutions and the SU(2) orientation of the vector q. 

Because the moduli space metric for the Dancer data is already known, the methods of 
Sec. II can be applied to obtain the metric for the full 40-dimensional moduli space. However, 
the result would be rather unwieldy for exploring the nature of the cloud dynamics. We 
will therefore reduce the problem to a more manageable one by restricting ourselves to a 
considerably smaller, but gcodcsically complete, submanifold. 

A geodesically complete submanifold can be obtained by restricting to the maximal sub- 
space left invariant by some isometry of the full manifold. In particular, we will require 
that the solutions be axially symmetric about the z-axis. This means that the Nahm data 
must be invariant under the combination of a rotational transformation of the form given in 
Eq. (2.2) and an appropriately chosen gauge action. In each set of Dancer data two of the 
fj must then be equal (which is only possible if = or 1), which implies that the solution 
acts as a symmetric top in the SU(2) space. Furthermore, the symmetry axes of the two 
Dancer solutions must be aligned with each other and with q. More specifically, the fj" and 
the fj^ can differ only by an U(l) rotation. Making use of the redundant U(2) freedom noted 
above, we can take the U(l) rotation to be about the direction and fix the fj" and fj^ to 
be 



Although rotation of the relative phase is not an isometry, there is a Z2 symmetry that 
reverses its sign. We can require invariance under this symmetry as well, and set ^ = 0.^^ 

If we now set = ~ and write T4 = p -|- qr-^^ the Nahm matrices on the left and 
right intervals then become 



T 



Ri 



\9i{s)ru \9ir2{s), ]^gi{s)n + ZRh 



(4.7) 



Invariance under this Z2 symmetry could also be achieved by setting V' = — 7i'/2; we will not explore this 
possibility here. 
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with 



-Dl csc[Dl{s - sl)] , kl = 
-Dl cosech[DL{s - sl)] , = 1 

-DLCot[DLis- Sl)], kl = 
-Dl coth[DL{s - Sl)] , = 1 
Dr csc[Dr{sr - s)], Kr = 
Dr cosech.[DR{sR - s)] , hr^I 
Dr coi[DR{sR - s)] , = 

Dr coih[DR{sR - s)] , kr^I . 

We can further simphfy matters by requiring that the conserved charges from the unbro- 
ken U(l) xSU(4) xU(l) symmetry all vanish. One's first thought might be that the phases 
associated with these vanishing charges could be simply dropped from the Lagrangian. This 
is not so, because there are couplings between the angular velocities a;* of these phases and 
the six non-phase moduli {Dl, Dr, p, q, Zl, and Zr) that remain after our symmetry con- 
straints are imposed. If we denote the latter moduli by j/", the moduli space Lagrangian can 
be written as 



(4.8) 



1 



^MS = ^Cab rr + Ba^ + -E,, Uj'uj^ , (4.9) 

where the metric coefficients Cab, Bai, and Eij depend only on the y". By means of a Lcgcndre 
transformation we can convert this to an effective Lagrangian in which the dependence on 
the a;* is replaced by a dependence on the conserved charges 



If all of the Qj vanish, this effective Lagrangian reduces to 



^MS.eff 



Cab — BaiE,/B 



5j6 



(4.10) 



(4.11) 



As we did for the (1,[1],1) example, we will take advantage of the isometries of the moduli 
space and calculate the metric at the point V — I. We start our calculation by displaying 
the Nahm data. The T^{s) and T^{s), as well as T4, were given above. With V — I, 
A = i^V2^ where 

fp + q + C + R \ 









p-q-C + R 2B 

2B p + q-C-R 







with 



B = 
C = 





9i{so) - 9i{so) 
93(^0) - 9i{so) 







p-q+C-Rj 



(4.12) 



R = Zr-Z 



R ■ 



(4.13) 

Note that K has been written so that the Greek indices in Eq. (2.6) label 2x2 blocks; the 
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elements within each block are labeled by the indices r and s. 

Given this Nahm data, the calculation of LMs.cff can be organized as follows: 

1) Calculate the derivatives of the Nahm data with respect to the On the 
left and right intervals the only nonvanishing derivatives are those of the with respect 
to the corresponding D and Z, but A has nonzero derivatives with respect to all of the y". 
The calculation of these is somewhat involved, and so we relegate it to the Appendix. The 
explicit form of the results are actually not needed until the final step 6. 

2) Determine whether the tangent vectors obtained in step 1 require any 
additional gauge actions to put them into background gauge. It is easy to see that 
the derivatives of the on the left and right intervals obey Eq. (2.11). The pieces arising 
from the data at Sq require a bit more care. With V taken to be the identity, A is Hermitian. 
Equation (2.12) then implies that a compensating gauge action is only needed if 



— A 

dr 



7^0 (4.14) 



for any coordinate y". To see that this quantity always vanishes, first note that both A 
and dA/dy"' have the same block diagonal form as K, and that a nonvanishing commutator 
can only arise from the middle 2x2 block. Within this block, the matrices are all linear 
combinations of the identity and the Pauli matrices and p^. Any commutator must then 
be proportional to Py, and thus would not contribute after the trace over a was taken. 

3) Determine which Bja are nonvanishing. Because the tangent vectors for the 
do not require compensating gauge actions, Bj^ is given just by the first term on the last 
line of Eq. (2.20). This gives 



Bja = — 27ri Tr 



— A 



(4.15) 



where tj is the Hermitian generator corresponding to the jth phase. From the remarks of 
the previous paragraph, we see that we can choose the tj so that the only nonzero Bja come 
from the generator that has a py in the middle 2x2 block and zeros elsewhere; we label this 

generator 

4) Show that the tangent vector corresponding to the U(4) action generated 
by t2 does not need a compensating gauge action. Referring to Eq. (2.12), we see 
that this is equivalent to showing that 

= (^^2^),,,,, (4.16) 

for all values of r and s. It is easy to verify that this follows from the symmetric block 
diagonal form of A. 

5) Calculate E^-i- Because the Bja vanish if j ^ 2, we only need this one element 
of the matrix E"^ . Using the fact that the ^2 tangent vector has no compensating gauge 
action, Eq. (2.20) gives 

E-^j = 27rTr {A{t2,tj}A) = 27rTr {K{t2,tj}) . (4.17) 



This cancellation is a consequence of the axial symmetry, because otherwise there is also a py contribution 
to K. 
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This vanishes unless j — 2, implying that 

6) Evaluate the Cab and the B20 and substitute the results into Eq. (4.11) to 
obtain IvMs.eff • The details of this are given in the Appendix. Instead of writing the result 
directly in terms of the y"', it is more convenient to express it in terms of the four eigenvalues 
of K, 

Ai = p + q + R + C , 
A2 = p- q- R + C , 



A+ = p-C+J4B^ + {q-R)\ 



and the variable 



We can then write 



A_ = p - C - JAB^ + {q- Rf , (4.19) 



iMs,.« = Mal + Mh zl + i)i + i)| + 1 g + li^^-^e^ , (4.21) 

where Idd is the function given in Eq. (A. 9). (Of course, when obtaining the equations of 
motion from this Lagrangian one must remember that the A^- and 6 are not independent 
variables.) 



C. The large- mass limit 

Considerable simplification can be achieved by working in the "large-mass limit" in which 
the massive monopole core radii, M^^ and are much less than all other relevant 

distance scales. There are four possible cases, depending on the values of kl and kr. We 
will examine the two with kl — kr- 



1. Hyperbolic solutions, kl = kr = 1 

Here we take fiL = MlDl/Att and ^r = MrDr/4tt both large, with Dl and Dr held 
fixed. In this limit the Dancer clouds have minimum size and Dl and Dr are the separations 
between the massive monopoles of the same species. Up to exponentially small corrections, 

B^O, C^~(Dl + Dr). (4.22) 



It must be kept in mind that this hmit involves a comparison between the monopole masses and the cloud 

sizes and massive monopole separations. While the masses are, of course, constant, the evolution of the 
other quantities may invalidate this limit at large times. This would happen, for example, in a geodesic 
motion that started with a large-mass k = Dancer solution, passed through the spherically symmetric 
point where the symmetry axes in Fig. 1 meet, and then moved out toward the large-mass k = 1 solutions. 
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Substituting these values, as well as the asymptotic values of I^^ and from Eq. (A. 11), 
into Eq. (4.21) gives 



^MS,eff 



27r(Zi - Z4) 



ip - q) {Zi - Z4) 



+ 



Ti{p - q) 



1 

2' 



\p-qf + [Z,-Z,f 
27r(Z2 - Z3) 



+ ^M. Zl + -M. Z| H- (p + q) (Z. - Z3) 



+ 



^2 ^ 2 " 3 

7r(p + q) 



where 



[(p+g)2-(Z2-Z3)2] 

^1 = 
Z2 = 

^3 = 
Z4 = 



.2 _ (^^ _ Z,f^ 

{p + qf + {Z2-Z^f 



(4.23) 
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Zl- — 
Zr + 
Zr- 



Dl 
2 ' 

Dr 

2 ' 

Dr 

2 



(4.24) 



Examination of Eq. (4.23) shows that the metric is the sum of two independent pieces, 
one involving Zi, Z4, and p — g, and one involving Z2, Z^, and p + q. Each of these describes 
a (1,[1],1) SU(4) system. [Indeed, this could have been foreseen by recalling the results of 
Ref. [11], where it was shown that the SU(6) solutions with two minimal Dancer clouds 
and all SU(2) orientations aligned were essentially superpositions of two independent SU(4) 
(1,[1],1) solutions.] The splitting of the metric here implies that the two nontrivial clouds 
are completely decoupled from each other. Hence, this limiting case does not shed light on 
the interactions between clouds, which is our primary interest in this paper. We therefore 
turn to the second limiting case. 



2. Trigonometric solutions, kl = kr = 

For these, we take /il = (sq - sl)-Dl = MlDl/At: and hr = {sr - So)Dr = MrDr/At: 
to be just less than the maximum allowed value, tt. In this regime the approximate radius 
of the Dancer cloud is 

a = — > M-^ . (4.25) 
2(7r-/x) ^ ' 

To leading order, then, we can write 

C ^-B ^{aL + aR) = a. (4.26) 
In addition, using Eq. (A. 10), we find, again to leading order, that 

I^j, dDl + I^j, ciL>2 ^ 47r ^ + 47r ^ = IGtt (d^/Kf + IGtt ~ad(t? , (4.27) 

dL 0.R ^ ^ 
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where (f) — tsirr^{^ ai/ clr)- 

We can take the center of mass, M^Zl + MjiZ^, to be at rest and define a reduced mass 
M. = MlMr/ [Ml + Mr). The efi^ective moduh space Lagrangian of Eq. (4.11) then reduces 
to ^ 

LMS.eff = Mie + ^Y.^ + + 47r i + IGtt a 0^ (4.28) 

Note that, except in the 0^ term, the Dancer cloud size parameters aL and only enter 
the Lagrangian through their sum a. This is a consequence of our having aligned the U(l) 
phases of the two Dancer clouds, as described in Sec. IV B. 

Finally, in the limit of large monopole mass wc can treat R = Z^ — Zr as being constant 
in time, and so drop the first term on the right-hand side of Eq. (4.28). For the sake of 
simplicity, we will set R — 0. In the large-mass limit in which we are working, this makes 
the system essentially spherically symmetric. It also sets 9 — —ta,n~^{2d/q). 



V. CLOUD DYNAMICS 

We now focus on the dynamics of the trigonometric solutions discussed at the end of the 
previous section. We work in the large-mass limit with R — 0. The eigenvalues Act of the 
matrix K are then 

Ai = p + q + d, 

A2 ^ p-q + d , 

A+ ^ p - d + \Jq^ -\- 452 , 

A_ = p-d- yjq^ + Aa^ , (5.1) 

As noted in Sec. II, these eigenvalues must all be positive. Applying this constraint to the 
smallest eigenvalue, A_, gives the inequality^'^ 

p-\q\>d>Q. (5.2) 



A. Asymptotic behavior 

The system is particularly easy to analyze at large times (either positive or negative). 
The eigenvalues are then all large, with p^q':^ d — aL-\- qr, and 

Ai pa X+Ri p + q, 

A2 ~ A_ pa p — g . (5.3) 

Substituting these into Eq. (4.28), and noting that the 6'^ term in the Lagrangian is sup- 
pressed, we see that the dynamics is well described by the Lagrangian 

i^,„ = ,(P±i)! + ,(|Lli)! + 43 + 43. (5.4) 

p+q p-q ttL ttR 



For a fixed static solution, q is naturally defined to be positive. However, when describing time-dependent 
solutions it is convenient to allow q to change sign when it goes through a zero. 
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This can be viewed as describing a system composed of four noninteracting spherical 
clouds: two "SU(4) clouds", with cloud parameters {p + q) and {p — q), and two Dancer 
clouds, with cloud parameters gl and gr. (We will refer to these cloud parameters as 
radii, but it should be kept in mind that the cloud structure does not allow a precise and 
unambiguous definition of its radius.) These evolve according to 



p±q = 



(5.5) 



where the Ci and ti are arbitrary constants. The total energy is divided into four separately 
conserved parts. 



E. 



71 



E, 



p + q 
{p - qf 



p-q 

El 
Er 



TT- 

p-q 
47r— = 

4.^ = 



ttC, 



7rC_ 



47rCi^. 



(5.6) 



Note that this asymptotic separation into noninteracting clouds did not require that 
(p + g)-cloud and the (p — g)-cloud be very different in size, but only that they both be much 
larger than the Dancer clouds. This can be understood by recalling the description of the 
corresponding static solutions in Ref. [11]. By analyzing the magnetic field in the regions 
between the cloud radii, it was found that the non-Abelian part of the effective magnetic 
charge, Qna, in each of the regions can be diagonalized, with^^ 



Q 



NA 



fdiag(0, 0,0,0), 
diag(0,-1.0. 1), 
diag(-l, -1,1,1) 
diag(-2, 0,1,1), 

Uiag(-2,0,0,2), 



r > p + g 

p — q^ r aR 
ur:^ r ul 



(5.7) 



In other words, the clouds act as if they have magnetic charges 



Qp+q 
Qp-q 
Qan 



diag(0, 1,0,-1), 

diag(l, 0,-1,0), 
diag(l, -1,0,0), 
diag(0, 0,1,-1). 



(5. 



These formulas imply that at very large times the clouds would be expanding at speeds greater than 
that of light. A more detailed analysis of cloud behavior [16] shows that at these times the moduli space 
approximation breaks down, and that instead the cloud expansion is best described as a wavefront moving 

at the speed of light. 

We have arbitrarily chosen gr > and an > ul- 
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(a) 



(b) 



FIG. 2: Two examples of cloud collisions. The horizontal axis represents time, and the vertical 
axis cloud size. The incoming p — q and p + q clouds, represented by the dotted and dashed lines 
respectively, collide with the Dancer cloud (solid line) and then expand to infinity. 

Thus, the (p + q)- and the (p — g)-clouds he in mutually commuting SU(2) subgroups of the 

unbroken SU(4), and so can only affect each other via interactions mediated by one or both 
of the Dancer clouds. When p±q ^ a, these interactions are negligible, in accordance with 
the discussion in Sec. HE. Similarly, the two Dancer clouds decouple from each other in 
this asymptotic regime, regardless of their relative sizes. 

B. Scattering 

We have a system of four clouds that are asymptotically noninteracting. The asymptotic 
solutions indicate that they are all contracting at large negative times, and expanding at 
large positive times. Their interactions at intermediate times can be viewed as a series of one 
or more scattering processes. These can be studied by starting with an initial configuration 
containing well-separated clouds and then, using numerical simulations, letting the system 
evolve under the equations of motion that follow from the Lagrangian of Eq. (4.28). 

We show two typical examples of this in Fig. 2. Both of these simulations were performed 
with the constant of motion J — a^cf) set equal to zero, so that the ratio of the Dancer cloud 
radii remains constant throughout. The evolution does not depend on this ratio, but only 
on the sum of the Dancer radii, a, which is shown by the solid line in these plots. There is 
some ambiguity in defining the size of the two SU(4) clouds [e.g., the differences between Ai, 
A+, and p + q are negligible at large times, but not necessarily when the SU(4) and Dancer 
clouds are comparable in size]. We have, somewhat arbitrarily, chosen to plot p-\- q (dotted 
line) and p — q (dashed line). 

These plots show several features, common to all of the examples that we have examined, 
that should be noted. First, the SU(4) clouds always remain larger than the Dancer clouds 
(in fact, larger than the sum of the Dancer radii), as should be expected from the bound in 
Eq. (5.2). In the asymptotic solutions, the SU(4) cloud radii have parabolic dependences on 
time, with a minimum radius of zero. In the actual interacting solutions, their behavior is 
rather similar, except that the vertex of the parabola is raised so that it occurs at or near 
the point when the SU(4) cloud radius is equal to a. (Given the ambiguity in defining the 
cloud radii, the distinction between exact coincidence of these values, as in Fig. 2a, or a 
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slight gap between them, as in Fig. 2b, is not meaningful ) In particular, the overlap (or 
near overlap) between the SU(4) clouds and the Dancer clouds is relatively brief, suggestive 
of a rather short and sharp interaction. 

Also, from examining simulations for a variety of initial conditions, we see that, just as 
in the asymptotic limit, the (p + q)- and (p — gf)-clouds do not appear to interact directly 
with each other. This suggests that we focus on the interaction of just one of these SU(4) 
clouds with the Dancer clouds. We can do this by choosing initial conditions such that the 
(p + ?)-cloud is very large (and therefore essentially not interacting with the rest of the 
system) at the time that the {p — q)- and Dancer clouds interact. In fact, we can simplify 
our analysis by taking the {p + g)-cloud to be at infinity; i.e., by taking the limit p oo, 
with p^ /p and 6 = p — q held fixed. In this limit Ai = 2p — d + d and A+ = 2p — S — d + 0{l/p) 
tend to infinity, while 



A2 = S + d, 
A_ = S — d . 



(5.9) 



If we drop the terms proportional to that decouple from everything else, and restrict our- 
selves to the J = case where the ratio of 0^/07? remains constant, the effective Lagrangian 
of Eq. (4.28) reduces to 



7r(5 + a)^ -K(5 — dy d 
2 [0 + a) 2 [0 — a) a 



(5.10) 



Keeping in mind that we expect 5 ^ a at large times, and noting that this purely kinetic 
Lagrangian is equal to the energy, we can write 



E 



7raa6 



52 - d' 52 
Es + Ea- 



d' 



+ 



A ~2 

Ana 



Tx5a 



naad 



(5.11) 



where we have defined SU(4) and Dancer cloud energies whose asymptotic values at large 

\t\ are 



52 
a 



(5.12) 



It follows that the trajectories at large negative times, when 5 ^ d, are of the form 



d{t) 



E, 



167r 



- (t - to) 



(5.13) 



The trajectories at large positive times are of the same form, except that the values of the 
various constants of motion are changed as a result of the interactions between the clouds. 

The form of Eq. (5.10) is strikingly similar to that of Eq. (3.13), with the Dancer cloud 
parameter d playing a similar role to the separation between the massive monopoles in 
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(a) (b) 

FIG. 3: A typical collision between the Dancer cloud and an SU(4) cloud. The horizontal axis 
represents time, and the vertical axis cloud size. The actual cloud trajectories are shown as solid 
lines. In (a) the dashed lines indicate the initial and final asymptotic trajectories of the Dancer 
cloud, while in (b) they indicate the asymptotic trajectories of the SU(4) cloud. 

the SU(4) (1,[1],1) solution, and the fixed monopole reduced mass /i replaced by the variable 
a~^. This seems surprising, since the previous case involved massive monopoles hitting an 
ellipsoidal cloud at two distinct points, while in the present case two nested spherical clouds 
are meeting each other at all points. Nevertheless, the similarity in the Lagrangians suggests 
that the interactions should be similar. In particular, the analysis of the SU(4) dynamics 
in Ref. [13] found that the interaction between the cloud and the massive monopoles was 
relatively brief, taking place over a distance of order This suggests similarly brief 

interactions in the present case, with the interaction largely restricted to the time when 
5 — a is itself of order a. We saw some indication of this, with all of the clouds present, in 
Fig. 2. We illustrate this more clearly in the two-cloud case in Fig. 3, where we show the 
transition from the initial asymptotic trajectories to the final ones. 

Equation (5.13) suggests that an arbitrary solution depends on four initial constants, ta, 
ts, Ea, and Es- It is clear that time-translation invariance can be used to eliminate one of 
these. In addition, the Lagrangian of Eq. (5.10) has some interesting scaling properties. The 
only effect of the rescalings 

d — i> a' = Aa , 
S ^ 6' = X6, 

t ^ t' = Kt (5.14) 

is to multiply the Lagrangian by an overall factor of A/k^. Hence, given any solution of the 
equations of motion, these rescalings will generate a two-parameter set of solutions. Thus, 
to study the full range of possible solutions we really only need to vary a single continuous 
parameter, which we choose to be Ea/Es. (Note that applying the constraint 5 > a m the 
asymptotic region implies that Ea/ Es < 4.) Also, since the rescaling cannot reverse the 
time ordering, we must consider separately the cases ta — ts > Q and — < 0. 

The range of possibilities is illustrated in Fig. 4. If ta — ts > 0, the collapsing SU(4) 
cloud collides with the Dancer cloud while the latter is also collapsing. Three examples of 
this are shown in Fig. 4a-c, with the value of Ea/Es increasing from one to the next. In all 
three cases the SU(4) cloud loses energy to the Dancer cloud. In the last case, where Ea/Es 
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(a) 



(b) 




is initially close to its maximum allowed value, the SU(4) cloud loses so much energy that 
the inequality Ea/Es < 4 is temporarily violated. Because the cloud radii both increase 
like Et^, there must be a second interaction in which the Dancer cloud overtakes the SU(4) 
cloud and transfers back enough energy that the inequality is satisfied at large times. The 
crossover from the behavior shown in Fig. 4b to that in Fig. 4c occurs when Ea/Es ~ 2. 

In the borderhne case, ta — ts = 0, the two clouds arrive at the origin simultaneously, as 
shown in Fig. 4d. In this case the asymptotic solution of Eq. (5.13) is exact for all times, 
and no energy is exchanged between the clouds. 
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AE^/Eg AE5/E5 




FIG. 5: Energy change between the SU(4) and Dancer clouds. In (a) and (c), both clouds are 
contracting at the time of collision, while in (b) and (d) a contracting SU(4) cloud collides with an 
expanding Dancer cloud. In all four plots the dotted blue curve shows the transfer observed in the 
simulations. In (a) and (b) the solid purple curve shows the prediction for an elastic collision, while 
in (c) and (d) it indicates the prediction for an inelastic collision with A = —^{civ'j- + AvsiVai — 4:v'^^). 

Finally, we come to the case where ta — ts < 0. Here, the collapsing SU(4) cloud only 
reaches the Dancer cloud after the latter has already reached its minimum size and begun to 
expand. If Ea/Es < 4 is sufficiently large, as in Fig. 4e, the Dancer cloud loses some energy 
to the SU(4) cloud, but continues to expand, although at a reduced speed. (This is then a 
time-reversed version of a solution with ta — ts > 0.) However, if Ea/Es < 4 is small enough, 
as in Fig. 4f, the collision can reverse the expansion of the Dancer cloud and have it shrink 
to zero radius a second time. The boundary between these two regimes is at Ea/Es ~ 2.6. 

While these plots are sufficient to provide a qualitative understanding of the interactions, 
it would be nice to have some more quantitative results as well. Let us first consider the 
energy transferred during the collision. The fact that the interaction between the clouds 
takes place over a relatively short time interval suggests a naive model that treats the 
interaction as an instantaneous elastic collision of two rigid shells, with kinetic energy and 
radial momentum {J2a MaVa) conserved. Because the Dancer cloud has four times the kinetic 
energy of the SU(4) cloud for the same value of the velocity [see Eq. (5.12)], we treat it as 
having four times the mass. It is then a straightforward matter to calculate the fractional 
energy transfer. The result is compared with the actual data from numerical simulations in 
Fig. 5. We see that the model captures the important features of the collisions. It accurately 
predicts that if the two shells collide while traveling in the same direction, the faster one 
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(a) (b) 



FIG. 6: The parameter p, which is a measure of cloud thickness, as a function of the energy ratio. 
The result for two collapsing clouds is shown in (a), and for a collapsing SU(4) cloud and expanding 
Dancer cloud in (b). 



will lose energy. It also agrees with the data in predicting that in a head-on collision the 
SU(4) cloud will lose almost all of its energy for large values of Ea/Es. Finally, it correctly 
asserts that for a head-on collision there is a critical value of the initial energy ratio below 
which the direction of energy transfer is reversed. 

This model works better than one might have hoped, but there is no mystery as to why 
the predicted and observed energy transfer disagree. First, the cloud trajectories are only 
approximate, and are altered by additional interaction terms that only become significant 
when the cloud radii are comparable. Second, the interactions are not truly instantaneous, 
but occur over a finite time interval as the clouds move through one another. Let us modify 
the statement of conservation of energy of the clouds by including an inelastic term A, 
defined in terms of the initial and final cloud velocities by 

1^1 + ]^vl = Kjjf + Kjlf + A . (5.15) 

Because all of the terms in the conservation of energy equation are quadratic in velocities, we 
looked for an expression for A that was quadratic in the initial velocities and that provided 
good agreement with the observed energy transfer. By trial and error, we found that taking 
A = — |(3f|j + 4t>5jt>aj — 4t>^j) provides excellent agreement with the results obtained from 
numerical simulations, as can be seen from the plots in Fig. 5. The exact dynamics that 
give rise to this this formula are still unclear to us. 

We argued previously that the interaction between the clouds is largely restricted to the 
time when 5 — a was itself of order a; this gives us a measure of the thickness of the clouds. 
To describe this more precisely let us define the beginning of the interaction to be the time 
when 20% of the total energy has been transferred from one cloud to the other, and the end 
of the interaction to be the time when 80% has been transferred. We also define 5q and clq 
to be the values of these variables at the beginning of the interaction and 

p=^^. (5.16) 

Figure 6 shows p as a function of energy ratio for two different regimes. The left plot 
shows p for an interaction in which a collapsing SU(4) cloud overtakes a collapsing Dancer 
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cloud. Because the cloud velocities are equal when Ea/Es = 4, approaching this value 
from below corresponds to decreasing the relative velocity of the clouds. This plot therefore 
shows that p decreases as the relative velocity is decreased. We see that the two clouds 
can approach quite close to one another before exchanging a significant amount of energy 
if they are moving slowly relative to one another. The right plot is for an interaction in 
which a collapsing SU(4) cloud collides with an expanding Dancer cloud. In this case the 
relative velocity increases with E^/Es- We see that as this increases the clouds begin to 
transfer their energy sooner, and hence at a greater separation. The maximum value of p is 
about 0.17 in the former case and 0.15 in the latter. The similarity of the values for these 
two collision scenarios would seem to indicate that the cloud thickness is a relatively small 
fraction of the cloud radius, approximately (0.15 — 0.20)ao- 

The behavior described by these two plots suggests that the clouds act as dissipative 
media, and that when the SU(4) cloud moves through the Dancer cloud, the energy loss 
increases with the relative velocity of the clouds. This explains why when both clouds 
are collapsing and their relative velocity is small, they can come very close together before 
significant energy is transferred. In the other situation, where the clouds collide head-on, 
the energy transfer begins very quickly because the relative velocity is large. This is also 
consistent with the behavior of the inelasticity in the collisions that we found previously. 

VI. SUMMARY AND CONCLUDING REMARKS 

In this paper we have used moduli space methods to investigate the properties of the 
massless magnetic monopoles that arise when a gauge theory is spontaneously broken to 
a non-Abelian subgroup. We have shown how the natural metric on the Nahm data for a 
class of SU(2M-|-2) solutions with 2M massive and M(2M — 1) massless monopoles can be 
obtained from the metric of a simpler class of SU(M+ 1) solutions. Using this approach, we 
have explicitly verified for the SU(4) (1,[1],1) case that the moduli spaces for the Nahm data 
and for the BPS solutions are isomorphic, thus lending further support to the conjecture 
that such an isomorphism holds in general. We then applied this method to the problem 
of obtaining the metric for the SU(6) (2, [2], [2], [2] ,2) solutions from the (2,[1]) SU(3) metric 
studied by Dancer. This gave us an effective Lagrangian for a class of axially symmetric 
solutions. This Lagrangian was then used to study the interactions of the clouds that are 
the semiclassical manifestation of the massless monopoles. 

By examining explicit spacetime solutions, it has been known for some time that the 
spacetime fields evaluated at the cloud radius are not qualitatively different from those 
at points slightly further from or closer to the origin. One might therefore expect the 
interactions between clouds to take place as if the shells of these clouds were diffuse. However, 
our simulations show instead that the clouds interact more like relatively thin, hard shells. 
In the collisions between an SU(4) cloud and the Dancer clouds the energy transfer takes 
place over a short interval before and after the cloud radii coincide, suggesting an effective 
cloud thickness that is roughly 15-20% of the cloud radius. 

Some intriguing open questions remain. It is known that in Type IIB string theory one 
can interpret Dl-branes stretched between D3-branes as the analogs of massive magnetic 
monopoles. This suggests that massless monopoles should, in some sense, correspond to 
Dl-branes of zero length connecting coincident D3-brancs. It would be desirable to clarify 
these ideas, and to see if they would help explain the properties of the clouds that we have 
found. One would also hke to understand better the role of massless monopoles in the 
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electric-magnetic duality of = 4 supersymmetric Yang-Mills theory, where they should be 
the duals of the "gluons" , the massless gauge particles of the unbroken subgroup. We hope 
that our results will help shed light on these questions. 
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APPENDIX: CALCULATION OF LMS,eff FOR THE SU(6) EXAMPLE 



In this appendix we present the details of the calculation of the moduli space effective 
Lagrangian, Eq. (4.11), for the cylindrically symmetric SU(6) solutions of Sec. IV B. 

We begin by calculating the Cab- Because none of the y" tangent vectors require a 
compensating gauge action, Eqs. (2.19) and (2.20) for the metric reduce to 



Cab^^TT dsTi^^ + Anl (isTr^^ + 27rTr I -!^^ + 



so 



dy"- dy^ 



Qya Qyb Qyb Qya 



The only nonvanishing contributions from the integral over the left interval are 



and 



'-DD 



J so [dDr 



(A.l) 



(A.2) 



(A.3) 



(The mixed integral I^jj is zero because its integrand vanishes point by point as a result of 
the trace.) Noting that 



dZr 



Su3 h , 



we see immediately that 



Equation (4.7) implies that 



I^z = 87r(so - Sl) = 2Ml . 



^DD = 27r 



SL 



ds 



so 



' dgi 
dD, 



dD, 



(A.4) 
(A.5) 

(A.6) 



Defining u — s — so and referring to Eq. (4.8), we see that for the axially symmetric solutions 
with the prime indicating differentiation with respect to u. Hence, 
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^DD 



27r /■A^i/(47r) 



2n r 



du 



27r /•JWi/(47r) 

/ du 



{2gl + gl) + 2«(25i5; + ^3^^) + u\2g'^ + 5;^) 



+ — {2ugl + li'^^^g) 



(A.8) 



[In the second equality we have used Eq. (4.3) and its cychc permutations.] This is now 
easily integrated to give 



(Ml, . , . cos/iL \ 
— [jiL - sm hl cos hl) (tan Hl- — — , 

Ml ^ . , 1 \ /; , \ f cosh 



{Hl — sinh hl cosh /iL)(tanh/iL — hl 



^//^sinh //L^ 



Kl = 0, 



(A.9) 



where /il = MlDl/{^tt). 

The integrals on the right interval can be evaluated in the same manner, and give the 
same result, except for the replacement of M^, ji^, and by M/j, /i^, and k^, respectively. 

We need some limiting values oi Idd- For k = and close to tt. 



ttM 



while for K = 1 and large 



2(7r - At)3 

_M 



1 + 



1 



n — fjL 



1 + 



(A.IO) 



(A.ll) 



To calculate the contribution from A — K^/^, we recall from Eq. (4.12) that K can be 
written in the block diagonal form 



K 



/Ai 0\ 

K 

V A2 / 



(A.12) 



where the 2x2 matrix K can be expanded in terms of Pauli matrices as 

k={p- C)h + 2Bp, + {R- q)p, . (A.13) 

This can be rewritten as 

K^U-^PU, (A.14) 
where U — ex.pi9py/2 with tan^ = 2B/{R — q) and P is a diagonal matrix with eigenvalues 

X± = {p-C)± ^4^2 + {R-qy . (A.15) 

The square root of K is also block diagonal, with the middle block being K^^'^ — U'^P^/'^U, 
whose derivatives are 



(A.16) 



To calculate the metric we need 
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+ ^(a,^)Tr {daP'/'[P'/\py])-\{dJ){d,9)Tr[P'/\p,f 



(A.17) 



Because both P and daP are diagonal, the middle two terms on the right-hand side both 
vanish. The remaining terms give 



(A. 18) 



Adding to this the contributions from the corner elements of K gives 



IV daA d,A = Tr daK'/'^d.K'/'' = ^ 



4A. 



(A.19) 



where the X^r are the four eigenvalues of K. 

Combining this with our previous results, we obtain. 

Cab dy'^dy' = 2Ml dZl + 2Mr dZl + 7^^ dDl + /^^ ^^^^ 

^ daXadbK ' ' ' 

+ TT 



(A.20) 



The next step is to calculate the B2a- These only get a contribution from the boundary 
term, and are given by 
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—2m Tr 
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(A.21) 



With the aid of Eq. (A. 16), this can be rewritten as 

52a = T^daeTx{[[P^'\py\,P^'^]py} 



pl/2 



Py 



(A.22) 



We also need £'22^- Referring to Eqs. (4.18) and (A. 15), we see that 

1 



^22 



47r(A+ + A_) 



(A.23) 



Using these last two results, we can calculate the correction term that converts Lms, 
Eq. (4.9), to -LMS.eff, Eq. (4.11). The terms quadratic in daO combine nicely, and we find 
that 
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Cab — -6a2 -£'22^-626 



dy'^dy" = 2Ml dZl + 2Mr dZl + I^j, dDl + I^j, dDl 

dy^'dy^ . (A.24) 
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